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^1 "... one thing seems certain. That is that the work of 

OO , Elie Cartan on connections, holonomy groups, and 

homogeneous spaces is the source of all that is 



o 



Tij" ■ interesting in contemporary differential geometry. " 






^^ , K. Nomizu (cited in Guggenheimer ^ 

o 

r-| ; 1 Introduction 

(— I \ In 1901 Darboux revisited [7] the problem of mechanics considered originally by Bertrand 

P] in 1857 (for a Russian translation of the 1901 Darboux article refer to ^). Later on 
the results of this investigation due to Darboux were included almost verbatim in the 
classical text of Whittaker on analytical mechanics [SO]. To this day Darboux' approach 
to the problem, which thus came to prominence with the classical Whittaker's book, 
^ ■ is widely used in the mathematical physics literature on integrable and superintegrable 

^ . Hamiltonian systems (see, for example, [SSllMlEI)- Thus, the idea used by Darboux 

in [7j to solve the Bertrand-Darboux problem was generalized to study Hamiltonian 
systems defined in Euclidean spaces of higher dimensions (see |2Hj and the relevant 
references therein). In this view, one can assert that the classical Bertrand-Darboux 
problem is "that special case that contains all germs of generality" (in David Hilbert's 
words). Accordingly, the aim of this review article is to revisit the classical problem 
(which was only partially solved by Darboux, for a complete solution, using Darboux' 
approach, see P]) once again and present a solution utilising the fundamental ideas 
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of Cartan's geometry (see |16| and the references therein). The approach is based on 
both recent IHEOlEIllSlOHllBlIIllIinillSland classical [3 ED results that incorporate 
Cartan's geometry, in particular, - the method of moving frames |5l llUl[TTll241l25J . - into 
the study of classical Hamiltonian systems, including the Bertrand-Darboux problem 
that concerns classical Hamiltonian systems with two degrees of freedom defined in the 
Euclidean plane. 

More specifically, we target to compare, using modern notations and language, the 
two approaches to the problem originally considered by Bertrand and Darboux, namely 
Darboux' approach and the approach via Cartan's geometry. This will enable us to 
conclude that Cartan's geometry is the most general setting for solving the Bertrand- 
Darboux problem and its generalizations. 

2 A brief historical review of the Bertrand-Darboux prob- 
lem 

An adequate presentation of the Bertrand-Darboux problem should begin with a brief 
review of Liouville's paper of 1846 JSj to which its origins can actually be traced. Recall 
that Liouville in JH] studied canonical Hamilton's equations governing the motion of a 
particle on a curved surface defined by a metric in isothermal coordinates under the 
infiuence of a potential energy depending on coordinates only. He observed that if 
the metric and the potential of the corresponding Hamiltonian function (total energy) 
admitted special separable forms 

rrf , \{pI+pI) + C{u) + D{v) 



in some system of coordinates {u,v), the Hamiltonian system defined by (|2.H) could 
be solved by quadratures. Here A{u),B{v),C{u) and D{v) are some arbitrary smooth 
functions. Throughout this paper the configuration space of the Hamiltonian systems 
in question is a two-dimensional Riemannian manifold {M,g) (the pseudo-Riemannian 
case can be treated similarly [201 122| ISHj ) , and the phase space is the cotangent bundle 
T*M. Note that the metric of the kinetic part of 1)2. 1|) assumes the following covariant 
form 

ds'^ = {A{u) + B{v)){du^ + dv'^). (2.2) 

The form (|2.1j) of the Hamiltonian function also implies additive separation of variables 
for the associated Hamilton-Jacobi equation. The metric ()2.2() and potential of the 
Hamiltonian (|2.1j) are said to be in the Liouville forms. The converse problem was 
considered in 1881 by Morera j23j who showed that if a Hamiltonian system with two 
degrees of freedom defined by a natural Hamiltonian 

H{q,p) = \9''PrPj + V{q), i,j = 1,2 (2.3) 

could be solved within the framework of the Hamilton-Jacobi theory of separation of 
variables, then its metric g with the components g^^ , i,j = 1,2 and the potential V{q) 



assumed the Liouville forms with respect to special separable coordinates {u,v). Here 
q = ((?i,(?2)) P = {pi,P2) are the standard physical position-momenta coordinates. We 
note that this conclusion does not exclude the possibility that the Hamiltonian system 
defined by (|2.,3j) separates in some other system of coordinates with respect to which 
the metric g and the potential V{q) are not in the Liouville forms. Morera has also 
demonstrated that in the Euclidean plane E^ separation of variables for the Hamiltnonian 
1)2. 1|) in the Liouville form occurs in Cartesian, polar, parabolic and elliptic-hyperbolic 
coordinates. The equivalence established by Morera appears to be incomplete, since it 
is not clear where the separable coordinates (u, v) are to be derived from. The answer 
to this question, at least partial, was given by Bertrand and Darboux [21 E]. Recall 
that Bertrand considered Hamiltonian systems given by the natural Hamiltonian I\'2.'A^ 
and admitting an additional first integral satisfying certain conditions. In this view he 
considered the governing equations describing the motion of a particle in the Euclidean 
plane ¥? under the infiuence of the potential force defined by a function V{q) dependent 
on the position coordinates q = {q^ ,(f')- Thus, he assumed that the Hamiltonian system 
defined by (|2.3|) admitted a first integral of motion of the form 

F{q,p)=K'^p^pj + B''pk + U{q), i,j,k = 1,2 (2.4) 

and then looked for the potential function V{q) in ()2.3|) satisfying this condition. Bertrand 
considers the case of -B^ = B^ = and then shows that the vanishing of the Poisson 
bracket = {H, F} yields two further conditions, namely the Killing tensor equation 

[9,K]=0, (2.5) 

and the compatibility condition 

d{kdV) = 0, (2.6) 

where [ , ] denote the Schouten bracket |27|, the (1, l)-tensor K is given hy K = Kg ^. 

Remark 2.1. In what follows we shall always assume B^ = B^ = in ()2.4|) . 

The Killing tensor equation 1)2. 5(1 puts in evidence that the functions K"^^ = K^^{q) 
in (|2.4j) are the components of a Killing tensor of valence two defined in the Euclidean 
plane ¥?. More specifically, solving the Killing tensor equation 1)2. 5|) with respect to the 
standard Cartesian coordinates of the Eulidean metric 

g = diQdi+d2Qd2 (2.7) 

yields the following general form 

K = (/?! + 2/3492 + /369i)9i 5i 

+(/33 - f3m - 13^2 - I3miq2)di 82 (2.8) 

+ (/?2 + 2/3591 +/36g?)92 0^2, 

where denotes the symmetric tensor product, di = ^, 82 = -^ and the arbitrary 
constants Pi, . . . ,I3q are the constants of integration. In local coordinates the Killing 



tensor equation ()2.5() is an overdetermined system of PDEs and the formula (|2.8() demon- 
strates that the vector space of its solutions /C^(E^) is six-dimensional: dim /C^(E^) = 6. 
Bertrand next turns his attention to the compatibility condition ()2.6|) . In terms of the 
Cartesian coordinates q = (^1,^2) it assumes the following form: 

-^ - -TTY 1 (/569192 + /34gi + /3592 - /53) + 

7^^^(/36gi - /^eg? + 2/3492 - 2/3591 + /3i - /32)+ (2-9) 

dV dV 

9gi 9^2 

The linear PDE of second order (|2.9() was later called the Bertrand- Darboux PDE. Solving 
(|2.9|) for V^(q) amounts to finding admissible potentials of the Hamiltonian systems 
defined by ()2.3() whose integrability is afforded by the existence of first integrals H2.4() 
which are quadratic in the momenta. Bertrand restricted his attention to the solutions 
having the form 

V[{q^-uf + {q2-vfl (2.10) 

thus recovering known results due to Euler and and Lagrange about the motion of a 
particle attracted by two fixed centres according to Newton's law. 

The main goal of the 1901 paper due to Darboux |2j was to show that the linear 
PDE ()2.9|) could be solved in "full generality", thus yielding a general potential V of 
the Hamiltonian (|2.3|) that was compatible with the first integral (|2.4jl . Before solving 
the PDE (|2.9() . Darboux ingeniously observes that it can be simplified without loss of 
generality. Indeed, by rotating and translating the axes one can simplify the general 
form of 1)2. 8|) . thus bringing it to a certain canonical form. Note that this coordinate 
transformation has no effect on the dynamics of the system. In addition, Darboux 
excludes the case when [3q = 0, taking /Jg = 1/2 and also assumes that the Killing tensor 
after rotating and translating the axes assumes the following canonical form: 

K=(l3i + ]^qi\ diQdi- ^qmdi 082+ (f32 + ^qf] 82 82. (2.11) 

Note, that the coordinates 91,92 and the parameters /3i,/?2 in 1)2. 11|) are not the same 
as the corresponding coordinates and parameters in (|2.8j) . They are connected via the 
action of the isometry group (see ()4.31|1 below) that consists of rotations and translations 
of the underlying space E^. 



Remark 2.2. The canonical form of the Killing tensor (|2.1H) can be simplified further. 
Indeed, one can introduce the following Killing tensor, given hy K = K — (329'- 

K=(p,-P2 + ^qi) 81Q81- ^qiq28i Q 82 + \ql82 82, (2.12) 



where K is given by (|2.11j) and g denotes the contravariant metric of the Euclidean plane 
E^ given in terms of the Cartesian coordinates q = (91, 92) (|2.7|) . Indeed, in terms of the 



Hamiltonian H H2.3() and the first integral F H2.4() determined by the canonical Killing 
tensor (|2.11|) this is equivalent to introducing the first integral F which is quadratic in 
the momenta and given hy F = F — (32H. Now the difference /3i — /32 can be treated 
as the only parameter entering the formula (|2.12j) . Moreover, without loss of generality, 
one can always assume /3i — /52 > 0. 

In what follows, Darboux proceeds to find the unknown V by solving the PDE (|2.(i|) 
determined by the canonical Killing tensor H2.11() . using the method of characteristics. 
Eventually he introduces special coordinates (u, v) with respect to which the PDE (|2.(i|) 
is easily solved, yielding the general form of the potential V in terms of {u,v): 

C{u) + D{v) 

V[U,V) = 2 2 • (2-1^) 

u — V 

One immediately recognizes that the potential V H2.13() is in the Liouville form (|2.1|) . so is 
the corresponding metric. Hence, the Hamiltonian system can be solved via separation 
of variables in the corresponding Hamilton-Jacobi equation. Moreover, the separable 
coordinates (n, v) are the parameters of confocal conies related to the original Cartesian 
coordinates q = (^1,^2) as follows: 

qi = A; cosh n cos f, §2 = ksmhusvuv, (2-14) 

where I3i — j32 = 2/fe^, /3i,/?2 as in 1)2. 12() . More specifically, the parameters u and v 
represent two families of confocal ellipses and hyperbolas intersecting each other orthog- 
onally. In this view these orthogonal coordinates are called elliptic-hyperbolic. In view 
of the fact that they afford separation of variables in the corresponding Hamilton-Jacobi 
equation of the Hamiltonian system ()2.3() . we call them orthogonally-separable coordi- 
nates. In conclusion to this brief review of the results due to Darboux, we note that he 
essentially solved two problems: 

(1) First, under certain assumptions he employed an appropriate transformation that 
consists of rotations and translations of the space ¥? to reduce a given Killing 



tensor in the form (|2.8|) to its canonical form (|2.11j) . Note this transformation 
of E^ iduces the corresponding transformation in the vector space of Killing two- 
tensors /C2(e2) (see (I02|l below). 

(2) Second, he diagonalized the resulting canonical Killing tensor, thus producing the 
corresponding orthogonally-separable coordinates (|2.14|) generated by two families 
of confocal conies. 

Remark 2.3. The metric g and Killing two-tensor K that determine the Hamiltonian 
H (|2.3|) and first integral F (|2.4j) in terms of the orthogonally-separable coordinates 
(n, v) are given by the following diagonal forms: 

^ (du Qdu + d^Q d^) (2.15) 



A;^(cosh u — cos^ v 



and 

K = -—- -2 5-— (cos^ vdu Qdu + cosh^ ud-a d^) (2.16) 

fc^(cosn u — cos^ V) 

respectively, where 9„ = ^, 9^, = ^. 

In what follows, we shall describe the above problems solved by Darboux from the 
viewpoint of Cartan's geometry, as well as consider the cases that he excluded from his 
considerations. The following example (see j21j . as well as the relevant references therein 
for more details) presents an illustration of how Darboux' method can be employed to 
find orthogonally-separable coordinates (n, v) for a Hamiltonian system with two-degrees 
of freedom admitting a first integral of the form (|2.4j) and thus solve it by quadratures 
as shown by Liouville jlSj . 

Example 2.1 (2nd Integrable Case of Yatsun). Consider a Hamiltonian system 
with two degrees of freedom defined in E^ by the following Hamiltonian: 

1 2A 

+^qf+qiq^)-2{ql + q^). 

It is known that the Hamiltonian system defined by (|2.17j) is completely integrable if 
g^ = 2A admitting in this case the following additional first integral independent of 
(|2.17|) . which is quadratic in the momenta: 

F{Q,p) = ((l2 + l]pl-{2qi-l)q2PiP2 + {qi-l)qipl-3qt- 
2qlqi+ql + 6qf + 2q,q^-3qf. 
Observe that the Killing tensor K determined by (|2.18|) is given by 

K = (^l + qrjdiQdi+(^^q2-qiq2JdiQd2 ^^ 19) 

+ {-qi + ql)d2Qd2. 

We compare the general formula (|2.8j) with the formula (|2.19j) . and conclude that in this 
case /?! = J, P2 = Ps = P4 = 0, /Js = — ^ and /Je = 1- Thus, in this case the PDE 1)2. 9|) 
implied by the compatibility condition (|2.H|) is given by the following formula: 

d^V d^V\ , 1 ^ 5^^ . 2 2 3, 

(91^2 - 7;q2) + a a yi2 -qi + qi + 7)+ 



(2.20) 



dqj dq^J'^'^^ 2^^' ' dqidq2'^' .!■..' 4. 
dV dV , 1, 

dqi dq2 2 

We perform next the change of the variables qi ^ z + ^ ^2 — > y; to obtain: 

^2/ + 7nr(y -^ +l) + 1^3y-— 3z = 0. (2.21) 



g^2 Qy2 J dzdy dz dy 



Consider now the differential equation of the characteristics of (|2.21() : 

zyidy"^ - dz'^) + (z^ - y"^ - l)dzdy = 0. (2.22) 

Introduce again the new variables u := z^ and v := y^ to transform (|2.22|) into the 

following ODE: 

dv\ dv , , , , 

-r] u-v + —{u-v-l) = {), (2.23) 

du J du 

which is a Clairaut's equation, having a general solution of the form 

(m + l)(mz^ -y^) -m = (2.24) 

in the original variables z,y. Re- write ()2.24() in terms of a new parameter a to get 

z^ y2 

1 = 1. 

a^ o^ — 1 

Observe that the characteristic curves of the PDE are two families of confocal conies. 
Therefore taking the parameters of the confocal hyperbolas and ellipses as coordinates, 
we have 

z = ab, y=[(a2-l)(l-62)]i/2^ 

or 

z = cosh u cos V, y = sinh u sin v. 

Our next step is to write the PDE (|2.2H) in terms of new variables a and b to obtain: 

n2 2.d^y .dV dV 

which has a general solution 

C{a) + D{b) 

a^-b^ 

in the Liouville form (|2.1|) . Finally, we transform back to the original coordinates qi and 
q2 to find 

gi = i + coshncosu, ^2 25) 

52 = sinh u sin u. 

Therefore the separable coordinates are of the shifted elliptic-hyperbolic type. Using these 
new coordinates, one can solve the corresponding Hamilton- Jacobi equation by separa- 
tion of variables and ultimately solve the Hamiltonian system in question by quadratures 
(for more details see HIj). 



3 Geometry of the Killing two-tensors 

As is well-known, the geometry of the aforementioned four coordinate systems that afford 
separation of variables in the Hamilton- Jacobi equation of a Hamiltonian system defined 
by (|2..S)) in E^ and admitting a first integral of motion (|2.4)) is induced by the geometry 
of the corresponding Killing two-tensor that determines H2.4() . More specifically, the four 
coordinate systems, namely Cartesian, parabolic, polar and elliptic- hyperbolic, are gen- 
erated by the eigenvectors (eigenvalues) of the corresponding Killing two-tensor. Thus, 
in each case the two families of confocal conies are the integral curves of the eigenvectors. 
In this setting we exclude from our consideration the trivial Killing tensors, which are 
the multiples of the metric g and as such do not generate orthogonal coordinate systems. 
Such Killing tensors determine the Hamiltonian (|2.H|) and its multiples. Each non-trivial 
Killing tensor has (almost everywhere) distinct and real eigenvalues. Another important 
feature (known already to Jacobi) of the four systems of coordinates is that the Carte- 
sian, parabolic and polar coordinate systems are degeneracies of the elliptic-hyperbolic 
coordinate system. Thus in the most general case a coordinate system generated by a 
non-trivial Killing tensor has two focii (i.e., the points where the eigenvalues coincide). 
This is the case (considered by Darboux [3|) of the elliptic-hyperbolic coordinate sys- 
tem. When the distance between the focii is zero, one gets the polar coordinates. Next, 
the case when one of the focii goes to infinity, corresponds to the parabolic coordinate 
system. Finally, when both focii are at infinity, the coordinate system is Cartesian. 
Figure 1 depicts the four coordinate systems along with the corresponding coordinate 
transformations to the Cartesian coordinates. 
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I: Cartesian coordianes 
qi = u, q2 = V 




H: Parabolic coordinates 

Qi = l/2{u^ -u^), q2 = uv 





IV: Elliptic-hyperbolic coordinates 

qi = A; cosh n cos w, 52 = A: sinh ti sin w 



III: Polar coordinates 

qi = ucosv, q2 = usinv 

Figure 1: Families of confocal conies 



We conclude therefore that Darboux in his celebrated 1901 paper considered 
only the most general case, namely orthogonal separation of variables in the case of 
the elliptic-hyporbolic coordinate system, excluding any degeneracies. When taking into 
account all of the cases, including the degeneracies, one is confronted with an equivalence 
type problem, that is for a given Hamiltonian system defined by (|2.3|) in E'^, admitting a 
first integral quadratic in the momenta, the first step is to determine what type of the 
four coordinate systems the corresponding Killing two-tensor generates. Once it is done, 
the next step is to transform the Killing tensor and its orthogonal coordinate system to 
a canonical form. As we have seen in Example 12. 171 on account of the existence of a non- 
trivial potential V in (|2.8j) . the system of coordinates that afford separation of variables 
can be not in a canonical form. Thus the elliptic-hyperbolic system of coordinates that 
we have derived in Example 12. 171 is shifted (translated along the (71-axis), that is affected 
by the action of the corresponding isometry group. Mapping a given element of an 
equivalent class to the corresponding canonical form is done with the aid of the moving 
frames map. Therefore we conclude that the most natural framework for solving the 
equivalence and canonical forms problems is Cartan's geometry. This observation leads 
us to believe that the most general solution to the classical Betrand-Darboux problem 
can be given in the framework of Cartan's geometry. This is the subject of considerations 
that follow. 



4 Cartan's geometry and Bertrand-Darboux' problem 

Let I{¥?) denote the Lie group of (orientation-preserving) isometries of E^. It is a semi- 
direct product of 50(2) (subgroup of orientation-preserving rotations) and T2 (subgroup 
of translations). Recall, the action I[¥?) O E^ is transitive, yielding a description of E'^ 
as a quotient I{E'^)/S0^{2), x G E^. The action /(E^) O iC^iE"^) (a general element of 
the six-dimentional vector space /C^(E^) is given by the generic formula (|2.8|) ) is not tran- 
sitive. Solving the equivalence and canonical forms problem is this case is equivalent to 



analysing the orbits of the vector space /C^(E^) under the action of /(E-^). The discussion 
presented in Section|31suggests that there are four types of orbits generated by non-trivial 
Killing tensors of the vector space /C^(E^). We shall exclude from further consideration 
the O-dimensional orbits generated by trivial elements of /C'^(E'^) (since they do not gen- 
erate coordinate systems and as such have no relevance in Hamiltonian mechanics). The 
orbit space /C'^(E'^)//(E^) has a rather complicated structure, since it contains 0-, 1-, 
2- and 3-dimensional orbits (see below). It has a local differential structure under the 
quotient topology. In addition, the projection vr : /C^(E^) -^ /C^(E^)//(E^) is a smooth 
fibration. The right action I(E?) O /C^(E^) is (locally) free. Furthermore, for any point 
X G /C^(E^) the corresponding orbit Ox G /C^(E^)//(E-^) is an immersed submanifold 
in /C2(e2). Thus, ^ = (/C2(E2),7r,/C2(E2)//(E2),/(E2)) is a principle /(E2)-bundle. For 
each orbit Ox passing through a non-trivial element x S /C^(E^) we can sonsider the 
bundle of orthonormal frames of eigenvectors (eigenforms) of the elements of Ox- The 
isometry group acts transitively on the bundle of orthonormal frames of Ox- Thus, we 
can employ the underlying ideas of Cartan's geometry, specifically - the moving frames 
m,ethod IIUI 1111 1241 I25j , - to solve the equivalence and canonical forms problem for the 
action /(E2) 0/C2(e2). 

Let X G K?{¥?) be a non-trivial covariant Killing two-tensor K. Recall |1] that in 
the rigid m,oving frame of orthonormal eigenforms E^,E'^ of K the components of the 
metric g of E^ and the Killing tensor K are given by the following formulas: 

gab = SabE'' Q E\ Kab = \a5abE'' & E\ a, 6 = 1, 2, (4.1) 

where 5ab is the Kroneker delta, A^, a = 1, 2 are the eigenvalues of K. The dual vectors 
Ei,E2 are the eigenvectors of K. At a given point x G E^ the two sets {E^,E'^^ and 
{El, E2} form non-coordinate bases of the cotangent and tangent spaces respectively. 

Remark 4.1. In the context of orbit analysis, chosing a rigid moving frame is equivalent 
to chosing an appropriate cross-section (see below). Hence, this technique will yield 
canonical forms of the orbits of the action I(¥?) O /C^(E^). 

Now we proceed to introduce the basic equations of Cartan's geometry. The equations 
for structure functions C^ab, a,b,c = 1,2 are given by 

[Ea,Eb] = C'^abEc or dE'' = ^C\,E'AE''- (4.2) 

We intorduce the connection coefficients F as follows: 

VE^Eh = Tab'E,, VeX = -^cdE\ (4.3) 

where V denotes the Levi-Civita connection associated with the metric g. 

Remark 4.2. The choice of the connection is not arbitrary. As is well-known from 
Riemannian geometry, given a connection V on a manifold M one can parallel propagate 
frames. For any path r between two points of M parallel transport along r defines a 
linear mapping L[t) between the tangent spaces of two points. This linear map is an 
isometry if the connection V is a Levi-Civita connectiction. Clearly, the linear map L{t) 
induced by a Levi-Civita connection V maps orthonormal frames to orthonormal frames. 
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The vanishing of the torsion tensor T'^^c is given by 

T\c = Tfec" - r.fe'^ - C\c = 0, (4.4) 

while the componets of the Riemann curvature tensor R"'hcd with respect to to the moving 
frame are defined as fohows: 

-K bed — -C/ci db + J- dfe -L ce — -f^di cb — >- cb >- de ~ '^ cd>- eb , l^.bj 

respectively. We now define a one-form valued matrix cj'^b called the connection one-form 
by 

co\ := Ta^E'. (4.6) 



Further, we can define 



l^ab := gadl^ b- 



On account of the above the connection one-forms, ojab are obviously skew-symmetric. 
They satisfy Cartan's structural equations, 

dE^ + Lu^'b AE'' = 0, (4.7) 

dw^ + w"cAw'=fe = e''6, (4.8) 

where we have introduced the curvature two-form 0"^ := {l/2)R°'bcdE'^ AE . Taking the 
exterior derivative of (|4.7|) and ()4.8|) yields the first 

e^b AE^ = (4.9) 

and second 

(ie^ + w%Ae"fe-e%Aw"fe = o (4.10) 

Bianchi identities, respectively. In addition, for a (0, 2) Killing tensor K we have the 
Killing tensor equation: 

Kiab-c) = 0, (4.11) 

where ; denotes the covariant derivative defined by 

Kab;c '■= EcKab - K^b^ca - Kad^cb ■ (4-12) 

Note that the equation (|4.1H) is the covariant version of the Killing tensor equation 
given in terms of the Schouten tensor for contravairant Killing tensors ()2.5() . We now 
adapt the equations of Cartan's geometry listed above to the study of the vector space 
/C2(E2). Define Tabc ■= Qcd^ab'^, Cabc ■= QadC^bc and Rabcd ■= gaeR^bcd- In the case 
of a two-dimensional (pseudo)-Riemannian manifold (M, g) there are only two linearly 
independent connection coefficients rii2, r2i2 and one component of the Riemann tensor 
i?i2i2- For convenience we write 7 := rii2 and 5 := T2i2- Since in our case {M,g) = E^, 
the condition of fiatness i?i2i2 = applied to (|4.5|) in the new notations yields: 

Ri2i2 = -EiS + E2-f--/^ -6^ = 0. (4.13) 
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The equations H4.2|) and 1)4 .7(1 transform accordingly: 

[Ei,E2] = --iEi-5E2, (4.14) 

dE^ =jE'^ AE^, dE^ = 5E'^ AE^. (4.15) 

In a similar way, the Killing tensor equation ()4.11() yields the following equations: 

E2X1 = 27(Ai - Ai), ^lAi = 26{X2 - Xi), E^Xi = E2X2 = 0, (4.16) 

where (|4.4|) has been used. The conditions for orthogonal integrability for Ei and E2 
follow from 1)4. 15() : 

E''AdE'' = 0, a =1,2. (4.17) 

Therefore by Frobenius' theorem there exist functions /, g and variables n, v, such that 

E^ = fdu, E^=gdv. (4.18) 

In addition, applying [£'1, £^2] to the eigenvalues Ai and A2 we arrive at the integrability 
conditions 

Ei-i = -37<5, E25 = 37(5. (4.19) 

It is easy to show (for more details see ^) that 

^ df I dg 

Moreover, f = g and 

f{u,v)=A{u) + B{v). (4.21) 

Now to obtain the expressions for the canonical forms of the four orbits it is natural to 
consider the following three isometrically invariant cases 

Case 1 ■J = 6 = 4^ Xi and A2 are constant 

Case 2 7 = 0, 5 / (7 / 0, 5 = 0) <^ Ai is constant (A2 is constant) (4.22) 

Case 3 7^ 7^ <^ Ai and A2 are not constant 

The first step is to solve equations (|4.1(ij) in terms of the variables u and v, which 
reduces to finding the corresponding expressions for the eigenvalues Ai and A2. Upon 
integration, substituting the formulas (|4.2fl|) and l|4.21|) into (|4.16j) and then transforming 
the result to its contravariant form, we obtain the following fundamental formula (for 
more details see jl]): 

K = iiKc + e29, 4,^2 GM, (4.23) 

where 

^c = ,, , ^ p, , {B{v)du Qdu- A{u)d, d,) (4.24) 

A[u) + B[v) 
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and 

9 = Af ^]n( ^ (^" Qdu + d^Q d,) (4.25) 

A[uj + i3(v) 

for some smooth functions A{u) and B{v), which are the same as in l\2.2\i . The for- 
mula (|4.23j) is the solution to the Killing tensor equation in the rigid moving frames 
of eigenvectors (eigenvalues) of the unknown K. Note it does not depend on the cur- 
vature of the underlying space. Importantly, (|4.2.S|) represents the canonical forms of 
the orbits of the action /C^(M)//(Af), where {M,g) is a two-dimensional Riemannian 
manifold of constant curvature, /C^(M) is the vector space of Killing two-tensors defined 
on {M,g) and I{M) is the corresponding isometry group. To specify the formula 1)4. 23(1 
to the Euclidean plane E^, we employ the flatness condition (|4.13|1 . and then find the 
corresponding A{u) and B{v) in each of the three cases (|4.22|) . The problem reduces to 
solving the corresponding differential equations determined by (|4.13|) . Each of Case 1 
and 2 leads to one solution and Case 3 yields two distinct solutions. Thus, as expected, 
in total we arrive at four formulas for Kc corresponding to four differnt types of orbits 
of the orbit space IC'^{E'^)/I{E'^) 01: 

Cartesian (C) : K^' = du Q du, 
Polar (P) : K^ = d^ Q 9„, 

PB v^du Qdu- v?d^ d^ (4.26) 



Parabohc (PB) : K 
Elliptic-hyperbohc (EH) : K 



^u cos^ vdu Q du + cosh udy d^ 
A;^(cosh u — cos^ v) 



Accordingly, the metric (j4.25|) that appears in the general formula for the canonical 
forms ()4.23|) , in each of the four cases is given by 



Cartesian (C) : g = duQ du + d^ Q d^, 

1 

u 



Polar (P) : g^ = du Q du + ^dy @ du, 



ParaboUc (PB) : ,PB , a„ fl + 8 O 8. _ (^ 

Elliptic-hyperbohc (EH) : g^^ = duQdu + duQdu 

fc^(cosh u — cos^ v) 

We next compare the metric (j2.7|) given in terms of the Cartesian coordinates {qi,q2) 
with the corresponding four metrics given by (J4.27J) in terms of the coordinates (n, v) and 
employ the method introduced recently in T5^ to derive in each case the transformations 
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from the coordinates {u,v) to the Cartesian coordinates {qi,q2)'- 
Cartesian (C) : qi = u, q2 = v, 

Polar (P) : qi = ucosv, 52 = usinv, 

(4.28) 
Parabohc (PB) : qi = \/2{v? - u^), 92 = uv-, 

Elhptic-hyperbohc (EH) : qi = /c cosh u cos f, 52 = /csinhiisinf . 

Finally, we use the formulas (|4.28|) to transform back to the Cartesian coordiantes ((?i, 92) 
the canonical Killing tensors given by ()4.26|) : 

Cartesian (C) : K^ = ^lO^i, 

Polar (P) : K^ = q^di 5i - qmdi ^2 + q\d2 82, 

Parabohc (PB) : K^^ = gs^i ^2 - 2^192 82, (4.29) 

Elliptic-hyperbohc (EH) : Kf" = (/c^ + q^)di Q di - qiq2di 82 

+qld2&d2. 

Now, we can describe the canonical forms of the orbits of every type. For example, a 
general formula of the canonical forms of the orbits whose elements generate the elliptic- 
hyperbolic systems of coordiantes can be obtained by substituting into the formula 1)4. 23() 
the cooresponding expression for the canonical Killing tensor K^ given by (|4.29|) and 
the formula (|2.7j) for g. Therefore, in terms of the Cartesian coordinates (91,^2) the 
canonical forms of the elliptic-hyperbolic orbits are given by the following family of 
Killing tensors: 

K^^ = hKf" + £29, hJ2eR, (4.30) 

where K^ and g are given by (|4.29j) and (|2.7j) respectively. We immediately see that, 
for example, the Killing tensor 1)2. 12() belongs to the family ()4.3U|) for ^1 = 1, £2 = 0, 
Pi — h = ^^ and as such represents a canonical form of the corresponding elliptic- 
hyperbolic orbit (see also 



Remark 4.3. We note that as far as applications of these results to the study of Hamil- 
tonian systems are concerned, it suffices to study only the Kc part of the formula ()4.23|) . 
Indeed, if the Hamiltonian system defined by the Hamiltonian H (|2.3j) admits a first in- 
tegral F 1)2. 4|) . then it also admits a first integral of the form H.1H + i2F, £1,^2 S M and 
vice versa. 

We have solved the canonical forms problem for the orbit space /C^(E^)//(E^). More 
specifically, for each type of orbits we have derived the most general form of the canonical 
forms. However, as we have seen in Example 12.11 in applications, such as the Betrand- 
Darboux problem, non-trivial Killing tensors can appear not in their canonical forms 
(unless, of course, F = in 1)2. 3() ). Thus, one has to be able to solve the equivalence 
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forms problem, namely for a given non-trivial Killing tensor (which may not be in its 
canonical form!), find a way to determine invariantly what type of orbit it belongs to. 
When it is done, one has to develop a systematic way of transforming a given non-trivial 
Killing tensor to its respective canonical form, as specified by the formula (|4.23j) . To 
achieve these goals, we shall employ another version of the moving frames method that 
has been introduced and developed recently ^1 ^2 El]- In our case, the idea is to 
identify the bundle of frames of eigenvectors of non-trivial Killing two-tensors with the 
group action /(E^) O /C^(E^) (since the group acts transitively on the bundle of frames) 
and then work in the group. We begin by determining the action of the Lie group I(E^) 
in the vector space /C^(E^) given by a general element ()2.8() . This action is induced by 
the action /(E^) O E^ given by 

gi = qi cosps - 92 sinpa + pi, 

(4.31) 
92 = qi smps + q2 cospa + p2, 

where pi,P2 and ps are the parameters of the isometry group /(E^). Using 1)4. 31(1 and 
the standard tensor transformation laws applied to the general Killing tensor (|2.8j) . we 
arrive at the following formulas: 

01 = Pi cos^ p3 - 2/33 cos p3 sin ps + P2 sin^ ps - 2p2 P^ cos ps - 2p2 (3^ sin p^ 

02 = Pi sin^ p3 - 2/33 cos p3 sin ps + P2 cos^ ps - 2pi P5 cos ps + 2pi P4, sin ps 

+P6pI 

03 = iPi - P2) sinps cosp3 + /33(cos2p3 - sin2p3) + (piP^ + P2P5) cosp3 (4.32) 

+{piP5 - V2Pi) sinp3 - P6P1P2, 

04 = Pi cos P3 + P^ sin P3 - /?6P2 , 

05 = Pb COS p3 - Pi sin p3 - PePi, 

06 = P6- 

The formulas (|4.32)) were first derived in [SI and then independently rediscovered in |21j . 
Next, using the standard techniques of the Lie group theory, we derive the infinitesimal 
action of the group /(E^) in the vector space /C^(E^) given by the following generators 
of its Lie algebra [^ : 

f d d \ , „ ^ , d ^ d „ d 



It follows from the formulas (|4.33)) that the dimension of the orbits, which is determined 
at each point by the number of linearly independent vector fields ()4.33|) . varies from 
to 3. For example, when /3i = P2 and P3 = Pa = P5 = Pq = the dimension of the orbits 
is obviously 0. 
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To solve the equivalence problem completely, we follow the approach introduced in 
[H] (see also J19j). First, we choose a cross-section K that intersects 3-dimensional orbits 
transversally where the group action is free and regular: 

/C = {/33 = A = /35 = 0}. (4.34) 

Next, we observe that chosing a cross-section in this setting is equivalent to chosing a rigid 
moving frame in the previous considerations, where we employed a classical version of 
the moving frames method. Moreover, the intersection points are precisely the canonical 
forms whose coordinates are given by invariants of the group action /(E^) O /C^(E^), 
that is the functions of the parameters /3i, . . . ,Pq that remain unchanged under the group 
action (see |21j for more details). Next, we derive the moving frame map 7 : /C^(E^) -^ 
/(E^) for the normalization equations corresponding to the cross-section H4.34() : 

/33 = /34 = /35 = 0. (4.35) 

Indeed, solving (|4.35|) for the group parameters pi,P2 and ps, we get [Hj 

/35COSP3 -/34sinp3 



Pi 



Pe 



/34COSJ?3-F/?5Sinp3 

P2 = ^ , (4.36) 

1 2(/33/36 + /34/35) 

P3 = - arctan ■ 



2 (36iPi-P2)-Pi + Pl 

Note, the moving frame map 7 : /C^(E^) -^ /(E^) given by (|4.3Hj) maps elements of the 
corresponding orbits to their respective canonical forms. We first observe from 1)4.32(1 
that 

Ai = /36 (4.37) 

is an /(E^)-invariant of the vector space /C^(E^). Substituting the expressions for pi, 
P2 and p3 into the first two formulas in ()4.32|) . we arrive at two additional fundamental 
I (K"^) -invariants: 

A2 = I36il3i+f32)-l3l-l3l 

(4.38) 

A3 = mPl-P2)-Pi+Plf + KP6P3 + P4P5?. 

Recall that the fundamental invariants Ai and A3 were derived first in [25, and then 
rediscovered more systematically, using the method of infinitesimal generators, in J2J 
(see also for other methods). The fundamental invariants Ai and A3 can be used to 
distinquish between the orbits |^[^E]. The classification is as follows. 

Cartesian (C) : Ai = A3 = 0, 

Polar (P) : Ai / A3 = 0, 

Parabohc (PB) : Ai = A3 / 0, 

Elliptic-hyperbohc (EH) : Ai / 0, A3 / 0. 
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(4.39) 



The fundamental invariants Ai and A3 in the case when both of them are not zero, 
have an interesting geometric meaning. The distance k between the focii of an elhptic- 
hyperbohc system of coordinates generated by a non-trivial Killing tensor, is clearly an 
/(E^) -invariant. Hence /c is a function of the fundamenal invariants found above. The 
corresponding formula was found in [2Tj to be 

/AT 



k = ^. (4.40) 

Given a non-trivial Killing tensor K = K?{¥?), we compare it with the general form 
1)2. 8|) , obtain the values for the corresponding parameters /3i , . . . , /Jg and then compute 
the funamental I(E^)-invariants Ai and A3. The problem of the determination of the 
orbit type for K, using the classification (|4.l-{9|) . is then straightforward. If K generate 
an elliptic-hyperbolic coordinate system, then we compute the corresponding moving 
frame map (|4.,'ffi|) . find the group parameters pi, p2, Ps, and then subsitutte them into 
the formulas for the group action (|4.31|) . The resulting transformation maps K to its 
canonical form. 

Dealing with the remaining three cases, namely Cartesian, polar and parabolic, is 
a delicate matter. The group action I{¥?) O /C^(E^) is not regular globally. Thus, for 
example, the vanishing A3 = indicates that the dimension of the orbits drops from 
three to two. If, in addition, Ai = the dimension of the orbits drops further to one. 
On the other hand, comparing the expressions for the generators of the Lie algebra 
of the isometry group /(E^) (|4..S,3j) with the formulas for the canonical forms (|4.29|1 . 
we conclude that the Cartesian orbits are 1-dimensional, polar - 2-dimensional and the 
orbits of the parabolic and elliptic-hyperbolic types are 3-dimensional. Clearly, in order 
to derive moving frame maps for the 1-dimensional and 2-dimensional and 3-dimensional 
orbits, other than the elliptic-hyperbolic type, one has to choose different cross-sections 
(for more details, see [201 El)- 

Finally, we will solve the same problem considered in Example 12. II using the results 
obtained in this section. 

Example 4.1 (2nd Integrable Case of Yatsun). This time we will not find the 
solution by solving the Betrand-Darboux PDE (|2.9j) for the potential V. Instead we will 
deal with the Killing tensor ()2.19() only. The material presented in Section ^ suggests 
that once the Killing tensor of (|2.4|) is available, solving the Bertrand-Darboux PDE 
1)2. 9|1 for V is redundant. Indeed, now we compare the Killing tensor 1)2. 19(1 with the 
general form (|2.8j) to obtain 

Pi = l, 132 = 0, (33 = 0, l3i = -\, /95 = 0, /36 = 1. 

Substituting this data into the formulas for Ai 1)4. 37() and A3 ()4.38|) . we obtain 

Ai = 1 / 0, A3 = ^ / 0, 
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which immediately shows that the Kilhng tensor H2.19() generates the elhptic-hyperbohc 
coordinates. Moreover, by (|4.4()j) the distance between the focii k = 1. Next, we compute 
the moving frames map H4.36() : 

1 

Pl = -^, P2=P3 = 0. 

Taking into account the formulas ()4.31|) {qi = Q^i + 1) and H4.28() (for (EH)), we conclude 
that the transformation to separable coordinates is given by (|2.25j) as expected. Thus, we 
have solved the same problem as in Example 12.11 without solving the Bertrand-Darhoux 
PDE \2.y\) ! Moreover, the procedure based on the moving frame method is purely 
algorithmical and so it can be implemented in a computer algebra package jl^j . 

5 Conclusions 

The analysis presented above (compare Examples 12.11 and 14. 1|) demonstrates that the 
Bertrand-Darboux problem (and its generalizations |14j ) can be solved within the frame- 
work of Cartan's geometry, in particular - its most important asset - the moving frames 
method. Moreover, the approach, based on the method of moving frames, is algorith- 
mical, independent of the curvature or signature of the underlying pseudo-Riemannian 
manifold {M,g), easily adaptable to different geometric settings. 
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